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Abstract 

The near-wall behavior of turbulence is re-examined in a way different from that pro- 
posed by Hanjalic and Launder^ and followers^ 3 ! ’M'l 5 !. It is shown that at a certain 
distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the 
smallest eddies in turbulence). All the important wall parameters, such as friction velocity, 
viscous length scale, and mean strain rate at the wall, are characterized by Kolmogorov mi- 
croscales. According to this Kolmogorov behavior of near-wall turbulence, the turbulence 
quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the 
large eddies become “ Kolmogorov ” eddies, can be estimated by using both direct numerical 
simulation (DNS) data and asymptotic analysis of near- wall turbulence. This information 
will provide useful boundary conditions for the turbulent transport equations. As an ex- 


ample, the concept is incorporated in the standard k-e model which is then applied to 
channel and boundary layer flows. Using appropriate boundary conditions (based on Kol- 
mogorov behavior of near- wall turbulence), there is no need for any wall- modification to 
the k-e equations (including model constants). Results compare very well with the DNS 
and experimental data. 

1. Kolmogorov behavior of near-wall turbulence 


It is well known that at a sufficiently high Reynolds number, R er , 
relations for the near-wall turbulence in a channel flow: 

there are following 

+ dU+ i 

—UV + + -j—r- — 1 

dy+ 

(1) 

In the inertial sublayer (y + >30), 


U + = j log y + + C 

k 

( 2 ) 

—uv + ~ 1 

( 3 ) 

In the viscous sublayer ( y + < 5), 


U + =y + 
—uv + ~ 0 

( 4 ) 

( 5 ) 


*Work funded under NASA Cooperative Agreement 3-233. 
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where 



uv^ — 


uv 



+ _ U rV 

y = — 


R e 


u T h 


v v 

and k = 0.4 is von Karman’s constant, u T is the friction velocity, h is the channel half 
width, y is the distance from the wall, U and uv are the mean velocity and turbulent shear 
stress respectively. The Kolmogorov length, velocity and time microscales of turbulence 
are 


t) = ( — )«, v = (ue)*, r = (-) 1 * (6) 

£ £ 

where v is the kinematic viscosity, and e is the dissipation rate of turbulent kinetic energy. 
Equations (l)-(6) will be used intensively in analyzing the near-wall behavior of turbulence. 


Now let us show that the important wall parameters, such as the friction velocity 
u T , the viscous length scale v/u T , and the mean strain rate at the wall, are actually 
characterized by Kolmogorov microscales. We may use Eq.(6) to write, 



where e + = ■ According to the direct numerical simulation dataW’^’W, the peak 

value of £ + near the wall increases slowly from 0.165 to 0.251 with increasing Reynolds 
number (based on the momentum thickness), and tends to an asymptotic value (see Fig.l). 
These DNS data show that is of order one. If we take 0.251 as the asymptotic value 
of at high Reynolds numbers, the value of is about 1.413. This means that very 
near the wall 

u r « 1.413 v (7) 

because vr\/v is always one. Equation (7) can be also shown in the following way: 

» = ( 8 ) 

where we have estimated £ with u\jky , which is accurate in the inertial sublayer, and may 
still give the order of magnitude of e further toward (but not at) the wall. Eq.(8) shows 
that v is indeed of order u T near the wall where ky + is of order 1. 


The viscous length is defined as u ju T . If u T is estimated with v, we obtain 

v v v 

« — = j- = 7] 

U r V (i/ £ )« 


( 9 ) 


The mean strain rate at the wall is of order dU/dy. Using Eq.(4), we obtain 

dU _ u 2 T ~ v 2 _ (i/£)? _ 1 
dy v v i sr 


( 10 ) 
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Equations (7)-(10) show that near the wall all important physical parameters of turbulence 
axe characterized by Kolmogorov microscales. We shall see later that there is a narrow 
region where the distance from the wall is smaller than the Kolmogorov length scale. In 
that region, the turbulence can not be self- sustained because of the overwhelming viscous 
action, and the flow is dominated by the viscous stress. Therefore, it is not important that 
what is the behavior of turbulence within that region ( y < 77 ). The important thing is 
the behavior of turbulence above that region. It is important to realize that approaching 
the wall, the turbulence is going to be described by Kolmogorov microscales , which are 
finite, and non-zero quantities. We will refer to this asymptotic property as Kolmogorov 
behavior of near-wall turbulence. 

2. Near-wall asymptotic turbulence quantities 

To study the near-wall behavior of turbulence, one needs first to find where (or at 
what distance from the wall) the size of energetic large eddies will reduce to Kolmogorov 
microscale, and second to estimate the turbulent quantities, such as the turbulent kinetic 
energy, the dissipation rate, etc., at that position. 

Let y v indicate the position where the integral scale of turbulence t equals the Kol- 
mogorov length scale 77 . We examine the variations of these two scales as the wall is 
approached from the inertial sublayer. The size of energetic large eddies, or the integral 
scale of turbulence is of order y. In the inertial sublayer, it can be estimated as (see 
Tennekes & Lumleyf 9 ! ) 

i^ky ( 11 ) 

For the Kolmogorov length scale 77 , if we estimate e with u\jky (which is quite accurate 
in the inertial sublayer), we obtain 


V ~ ( 



( 12 ) 


Equations (11) and (12) show that both t and 77 decrease toward the wall, but at different 
rates. The integral scale i decreases much faster than 77 does. We may imagine that at 
some point near the wall, i will approach rj. We refer to this point as the limit point of 
turbulence, because passing beyond this point toward the wall, the turbulence is no longer 
important. Now if we use Eq.s (11) and ( 12 ), and let i equal 77 , we obtain 


(Jfe^)f =1 


V 


(13) 


Equation (13) indicates that the limit point of turbulence is 



= 2.5 


The numerical value of here should not be taken too seriously. It is an order of 
magnitude value, because we have assumed e — u\jky which is accurate only in the inertial 
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sublayer, and also this value can be changed by the point we pick for the viscous cutoff in 
the spectrum. It is reasonable to choose the value at the top of the viscous sublayer: 

vi ~ 6 (14) 

In the region where y < y v , the self-sustaining turbulence cannot exist. Here we define 
turbulence to be a self-sustaining chaotic motion of the fluid, dominated by inertia, so 
that the scale of the dissipation is at least somewhat smaller than the scale of the energy 
containing eddies. As viscous effects become stronger, and the two scales approach each 
other, a point comes at which the chaotic motion can no longer sustain itself, but must 
depend on energy imported from neighboring regions. The chaotic disturbances of the 
viscous sublayer are of this nature; by our definition they are not turbulent, although 
surely chaotic, since they depend on the truly turbulent motions in the buffer layer and 
beyond for their energy. 

The physical importance of Eq.(14) is that there exists a turbulence limit point y v 
which is not zero. To study the asymptotic behavior of turbulence near the wall is to find 
the behavior of turbulence at this turbulence limit point. 

Now let us estimate the turbulence quantities at the point y v . First, let us look at the 
dissipation rate at that point, e v . In the inertial sublayer, e is well represented by u\/ky 
or ( u*jv)(l/ky + ) which shows that the dissipation rate increases toward the wall. If we 
extrapolate this relation to the limit point y v , we obtain that e v is of order u*/u. On the 
other hand, if we use Eq.(7) at y v , 


", = 1^3 = (**,)* ( 15 ) 

we obtain 

e„ = 0.251— (16) 

v 

This means that according to the DNS data, e v is indeed of order u*jv. Eq.(16) gives the 
peak value of the dissipation rate near the wall. 

At the turbulence limit point, the characteristic velocity is v v , therefore the turbulent 
kinetic energy K at y v can naturally be estimated as 



u 


2 4.413 ; 


0.250t4 


(17) 


For the turbulent shear stress — uv, if Eq.(l) and (4) are used at y n , we obtain 


-(««), = 0 


(18) 
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Other shear stresses should also vanish at y v . Finally, the energy components u 2 , v 2 , and 
w 2 are not equal because of the anisotropy of turbulence near the wall. If we use the 
experimental dataM i n the inertial sublayer to estimate the anisotropy at y v , i.e. 


U 2 

= = a w 6.250 
v 2 

w 2 

= = b w 3.625 


and u 2 + v 2 + w 2 — 2K , we obtain 


v 2 

w 2 


a 


a b \ 

1 

a + 6 + 1 

6 

a + 6 + 1 


2tf 

2K 

2K 


(19) 


( 20 ) 


Using Eq.(17), we may estimate these energy components at y n \ 


(u 2 )^ w 0.287 i4 

(v 2 ), « 0.046^ (21) 

(w*), w 0.167u^ 

Equations (14), (16), (17), (18), and (21) represent the near-wall asymptotic behavior of 
turbulence. These equations provide useful information about the boundary conditions 
for turbulent transport equations such as K-e equations, Reynolds stress equations, etc. 
These ideas can be used at any level of turbulence modeling. In the next two sections, we 
illustrate how they can be very easily incorporated into the eddy viscosity K-e models. 


3. Eddy viscosity 


Eddy viscosity (for incompressible flows) is defined as 

-uiu] = vr(Uij + Uj,i) - ^Siju k u k ( 22 ) 

For a two-dimensional channel flow, it becomes 

_ dU 

—uv = vt~j- (23) 

dy 

Note that Equation (1) is valid for any finite y + (but not for y + — » oo at R er —■ ► co), 
therefore we may use Equation (1) to study the behavior of the eddy viscosity vt- If we 
approximate dU / dy with u T jky, and use Eqs. (1) and (23), we obtain 

v T = ( ky + -1) v (24) 
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This equation shows that the eddy viscosity is ku T y for a large . Now if we model vy as 


vr = CJ f — (25) 

£ 

and use experimental data to estimate K » 3.5*4 in the inertial sublayer, and e ~ ul /ky, 
we obtain 

wt = 12.25 C tl f ll ky + u (26) 

Comparing (24) and (26), we obtain 

f, 

C, 


= 1 - 


ky+ 


12.25 


0.082 


(27) 


Eq.(27) shows that the damping function ~ 1 when y + > 30, which means that the 
standard K - e eddy viscosity model is reasonable in the inertial sublayer, and a significant 
damping is needed only below < 30, which corresponds to the buffer layer. However, we 
do not expect that the form of (27) will be appropriate for the buffer layer region, because 
the assumptions made in these relations are restricted to a large y + . Unfortunately, we 
are not able to find similar analytical expressions like Eqs.(24), (26) in the buffer layer to 
form the damping function like Eq.(27). We must then depend on the experimental or 
DNS data to propose appropriate form for C M and f ^ to fit Eqs. (23) and (25). Based on 
the work of Yang and Shih f 10 ^, we propose 


Cfi = 0.09 

/^ = [1 — exp(aii2)t + + agi2|)] 2 

ai = -1.5 * 10 -4 a 3 = -1.0 * 10~ 9 a 5 = -5.0* 10 -1 ° ( 28 ) 

R k = 

V 

This is no more than a curve fitting based on the solution of DNS of channel flow with 
R er = 180. Fortunately, it also works quite well together with the K-e equations (described 
in the next section) for other cases (see section 5). We choose Rk in Eq.(28) instead of y + 
in order to avoid the unphysical behavior of near the flow separation or reattachment 
points where the friction velocity u T is zero. 

4. K-e equations 

We may use standard k-e model equations for the eddy viscosity model (25), and use 
Eq.(22) to calculate turbulent mean flows. The boundary conditions can be obtained from 
Eqs. (14), (16), and (17). The K-e equations for incompressible flows can be in general 
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modeled as (see Appendix 1 for the derivation of e equation): 


^ = [{« + -)K,i},i + nVijWtj + Vi,,) - c (29) 

Dt <?k 

— =[(* + + CihvTVi,i(v,,i + &(,<)£ 

+vv T U,,i k U, Jk (30) 

These equations axe used only for the flow field outside of the turbulence limit point y v , 
where K v is non-zero. Therefore, Eq.(30) will not have singularity problems and will not 
need any near-wall modifications like other K-e models doJ 11 !^ 12 ! 

As examples, we write the equations and the boundary conditions for both the two- 
dimensional channel and boundary layer flows. 

4.1 Channel flows 


For channel flows, all the quantities are normalized by the friction velocity u r , and 
the half width of the channel h , and v? is normalized by v . 


dU 

dt 

dK 

dt 

de 

dt 


d .. 1 ,, , dU 1 


— [— (1 + — )— ] + 

dy V R eT K ^ <T k > dy ] ^ R tr 


1 ,dU, 2 


d , 1 VT.de, . 

o- -5- 1 + — )«- + Cifi 
dy R er <? t dy } 


1 dU . 2 e nr* 2 , 1 

R ct UT ^’ K~ 2 h K + R2^ 


UT 


a 2 ir 

dy 2 


(31) 

(32) 

(33) 


where R er is a known parameter for channel flows. The boundary conditions for the mean 
velocity are straightforward: 

QTT 

y = 0: U = 0, y = 1: — =0 (34) 

dy 

The boundary conditions for K and e must be given at y v . In the normalized form, 
y = y+ /R er . Using Eq.(14), (16) and (17) we obtain 


y =Vv 

y =1 : 


g 

— 1 Er) = 0.251 Ren Krj 

tier 

de „ dK „ 
dy ~ °» dy ~ 


0.25 


( 35 ) 


In solving Eq.(31), the eddy viscosity vt is negUgible when y < y v . 
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4.2 Boundary layer flows 


For a turbulent boundary layer flow, all the quantities are normalized by the free 
stream velocity Uqq , and the length scale L. ut is normalized by v, and the pressure P is 
normalized by pU 


DU 

Dt 

DK 

Dt 

De 

Dt 


d f 1 n , .dU 


y-[j-(l + =)^l + 

dy Roo <?k oy 
d , 1 „ , VT.de 


dP 

dx 

1 ,dU \ 2 

- c 


dU, 


ir.lihV + - c 2 f 2 - + 


dy R 


<t £ dy 


dy' K 


K 


R 2 

±L oo 


Vt 


d 2 U‘ 

dy 2 


(36) 

(37) 

(38) 


Boundary conditions for the mean velocity are: 


y = 0 : U = 0, 


y = OO : 



(39) 


The boundary conditions for K and e are given at y v . In the normalized form, y = 
y + 1 {RooUt)- Using Eq.(14), (16) and (17) we obtain 


V =Vv = 


e v = 0.251i2oottr> K v = 0.25u 2 


RooUr 

y =oo : 

where the friction velocity is calculated from the solution of the mean velocity. 


de n dK n 
dy ~ °’ dy ~° 


(40) 


J_dU. 
Roo dy j * =0 


(41) 


In solving Eq.(36), the eddy viscosity ut is zero when y < y v . 

In practical apphcations, R>c t and R-c oo are large numbers, hence y^ is usually very 
small. Therefore, as an approximation we may let y v = 0, but e n and K v must be given 
by Eqs.(35) and (40) respectively. These equations have been applied to the calculations 
shown in the next section. 


5. Comparison of models 

To compare the present model with the DNS data and other models (e.g. Jones and 
Launder^ 11 ], and Chient 12 !), we have made calculations on two channel flows t®ld 7 l and two 
boundary layer flows^’l 15 !. In the present model, all the model constants are the same 
as used in the standard K-e modefl 13 !. Therefore the present model will also be suitable 
for flows far from the wall. The other two models used here for comparison do not have 
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this property. Results are shown in figures 2-5. In figure 2 and figure 3, three models 
are compared with two DNS data for channel flows: one with R eT = 180, the other with 
R er = 395. The profiles of mean velocity, Reynolds stress, turbulent kinetic energy and its 
dissipation rate are plotted in these figures. The present model is significantly better than 
the other two models. Figure 4 shows the similar comparison for a turbulent boundary with 
R e0 = 1410. The agreement between the present model and DNS data is excellent. Figure 
5 shows the results compared with Klebanofft 1 ^ and other boundary layer experiments. 
The sVin friction from DNS data^l is also shown in this figure. The results of present 
model are more consistent with the DNS data than the experiments. 

It is also worthwhile to emphasize that the present model equations with the standard 
model coefficients have the simplest form among all two-equation models. Hence, we expect 
that they will have less numerical stiffness in complex turbulent flows. 


Appendix 1. K-e equations 


The equation for the dissipation rate of turbulent kinetic energy e is proposed following 


Lumleyt 14 !: 


De .. 

m= l(v + 




(Al) 


where $ stands for the entire mechanism of the production and destruction of the dissipa- 
tion rate e. At the level of the K-e model, we assume that $ is a function of i/, vt, K , e^Uij, 
and Uijk* Because ^ is an invariant, it must be a function of the invariants that can be 
constructed from these quantities: Rt^T^i^Ui^ / e, and uutV i } jkU ijk / £ 2 > where R t is the 
turbulent Reynolds number ~~ . We now expand ^ in a Taylor series about these invariants 
and keep only the linear terms to obtain 


^ = ^0+^1 


e 


+ VVtU i t jkU 


K 


(A 2) 


where the coefficients i>i and ip 2 are in general function of Rt- Inserting (A2) into (Al), 
we obtain 

^ = [(« + + Cl hvrViAVij + 

m a. K 

e 2 

- C 2 f2JZ + VV T Ui,jkUi,jk 

where Ci, C *2 are model constants, and /i ,/2 are in general functions of Rt , and ip 2 — — 1. 
The forms of Cj , C 2 and f\ , /2 are chosen to be the same forms as used in the standard 
K-e model: 

Ci = 1 . 44 , C 2 = 1.92 

/i = l> Si ~ 1 -0.22 exp (-”) (^ 4 ) 

(T k = 1 , <Tt = 1.3 
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Figure 1. The asymptotic value of the dissipation rate at high R e $. 
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Figure 5. Comparison of models with, the experiments of 2-D boundary layer 
flows (Klebanoftf 13 ] and others). 
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